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el Abstract. The microcanonical entropy s(e,m) as a function of the energy e 

' and the magnetization m is computed analytically for the anisotropic quantum 

Heisenberg model with Curie- Weiss-type interactions. The result shows a number 
of interesting properties which are peculiar to long-range interacting systems, 
including nonequivalence of ensembles and partial equivalence. Furthermore, from 

C^ , the shape of the entropy it follows that the Curie- Weiss Heisenberg model is 

indistinguishable from the Curie- Weiss Ising model in canonical thermodynamics, 

C^ ' although their microcanonical thermodynamics in general differs. The possibility 

of experimentally realizing quantum spin models with long-range interactions in 
a microcanonical setting by means of cold dipolar gases in optical lattices is 
discussed. 



1. Introduction 



^ ■ Many general theorems in statistical physics, but also the majority of applications, 

f^ I are concerned with short-range interacting systems. Short-range can either mean an 

'/^ ■ algebraic decay, r~''+", of the interaction potential with distance r, spatial dimension 

d, and a negative constant a, or a finite range of interaction. For proving theorems in 
statistical physics, short-range interactions are a very handy assumption, as they allow 
1^ I one to partition a large volume into smaller subvolumes while neglecting the surface 

effects of the subvolumes in the thermodynamic limit. This trick is at the core of the 
proofs of the existence of thermodynamic potentials and their convexity properties in 
the thermodynamic limit, and also of the equivalence of statistical ensembles like the 
microcanonical, the canonical, and the grandcanonical one. 

But there are also physical factors which account to some extend for the focus 



5— ( ' on short-range interactions: statistical physics was originally invented by Boltzmann, 

Gibbs, and others for studying the behaviour of gases, liquids, and solids, and the 
fundamental interactions in such systems are typically of an electromagnetic kind. 
The presence of positive and negative charges causes screening effects which give rise 
to interactions that are effectively of short range. For almost all practical purposes, 
effective interactions of finite range approximate this situation excellently, and this 
accounts for the almost exclusive interest in short-range interactions in the community 
of condensed matter physicists. 
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The situation is different in astrophysics where gravitational interactions are 
of relevance. Since masses are non-negative, screening effects do not occur, and 
gravitational interactions retain their long-range character. As a consequence, it was 
in the astrophysical context that peculiarities of the statistical physics of long-range 
interacting systems first attracted attention. Negative heat capacities were shown 
to exist in bounded self-gravitating gas spheres by Lynden-Bell and Wood [1] in 
1968, and Thirring explained this observation by relating it to the nonequivalence 
of microcanonical and canonical ensembles [2]. Although equivalence of ensembles 
had been proven only for short-range interactions, it was tacitly assumed by most 
physicists to hold in general. Therefore it came as a surprise to many that equivalence 
does not necessarily hold for long-range systems, in particular in the presence of a 
discontinuous phase transition. 

Within the condensed matter community, interest in long-range interacting 
systems arose somewhat later, and physicists with an interest in fundamental issues 
of statistical physics studied simple long-range toy models in order to explore 
the terra incognita beyond traditional short-range statistical physics. In order to 
study nonequivalence of ensembles, it became necessary to compute thermodynamic 
functions within the microcanonical ensembles. This is in general more difficult 
than computations in the canonical ensemble, but a number of exact microcanonical 
solutions of long-range systems have been reported (see [3] for a review) . The models 
used in these studies — most prominently among them the Hamiltonian Mean-Field 
Model [4] — consist of classical spin variables, and in most cases the interactions are of 
Curie- Weiss- type, i.e. each particle interacts with every other at equal strength. Curie- 
Weiss-type interactions can be considered as a limiting case of extremely long-range 
interactions, decaying as r~'^^°' in the limit a /^ d. 

For quantum systems with long-range interactions, only very few microcanonical 
calculations have been reported in the literature. A notable example is the paper 
by Pflug [5] on gravitating fermions, where a negative specific heat is found for all 
negative values of the energy. In the present article, we consider a different class 
of quantum systems, namely long-range quantum spin systems. The main goal is 
to contribute towards the understanding of nonequivalence of ensembles in quantum 
spin systems. The model chosen for this study is the anisotropic quantum Heisenberg 
model with Curie- Weiss-type interactions as defined in section [2] This is a model 
of anisotropically interacting spin- 1/2 degrees of freedom, and the Curie- Weiss- type 
interactions render it possible to find an analytic solution for the microcanonical 
entropy in the thermodynamic limit. It is one of the main purposes of this article 
to show how such an analytic calculation can actually be done, and the presentation 
in section [3] is therefore somewhat technical. In section 21 the canonical Gibbs free 
energy is recovered from the microcanonical entropy by means of a Legendre-Fenchel 
transform. Although the canonical result was known previously, its recovery from 
the microcanonical entropy is instructive and useful for discussing nonequivalence of 
ensembles and related issues in the subsequent sections. 

In the remaining sections, physical implications of the results are discussed. In 
section [SJ nonequivalence of ensembles is discussed. In particular, it is shown that 
microcanonical and canonical ensembles are nonequivalent for a certain range of 
anisotropy parameters of the Curie- Weiss quantum Heisenberg model. A related 
observation is that, for the same range of anisotropy parameters, the canonical 
Gibbs free energy of the Curie- Weiss quantum Heisenberg model is identical to 
that of a (classical or quantum) Curie- Weiss Ising model. These two models are 
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therefore thermodynamically equivalent in the canonical ensemble, but not so in 
the microcanonical ensemble, as discussed in section [51 The physical relevance of 
a microcanonical treatment of quantum spin systems is discussed in section [T] It 
is argued that the conditions under which a microcanonical description is adequate 
can be realized in experiments with cold atoms in traps. Dipolar gases in optical 
lattices can be used to engineer a large variety of long-range interacting condensed 
matter Hamiltonians, including the anisotropic quantum Heisenberg model 6 , and the 
absence of a coupling to a heat bath renders the microcanonical ensemble appropriate 
for a statistical description. We conclude in section [5] with a summary and a few 
further remarks. 

Some of the results reported in this article, and in particular a discussion of their 
relevance for experiments with cold gases in optical lattices, have been previously 
announced (without the details of the calculations) in a Letter [7]. 

2. Curie- Weiss anisotropic quantum Heisenberg model 

The model consists of N spin-1/2 degrees of freedom, each of which is interacting 
with every other at equal strength (Curie- Weiss-type interactions). The corresponding 
Hilbert space V. = ((D^)®^ is the tensor product of N copies of the spin-1/2 Hilbert 
space (D^, and the Hamiltonian operator is given by 

N N 

^" --m^ ^^'"'^"^ + ^"''^''' + ^'"''"''^ -hT.^1 (1) 

k,l = l k=l 

The cr^ are operators on % and act like the a-component of the Pauli spin-1/2 operator 
on the fcth factor of the tensor product space H, and like identity operators I2 on all 
the other factors, 

cr^ = I2 ••• (g) I2 «) cr" (8)l2 «) • • • (8) I2, a e {1, 2, 3}. (2) 

fcth factor 
The resulting commutation relation is 

[a,",af] =2i4,iea^W> a, /3,7 e {1, 2, 3}, (3) 

where 5 denotes Kronecker's symbol and e is the Levi-Civita symbol. The parameter 
h in the Hamiltonian is the strength of an external magnetic field orientated along 
the 3-axis, and the constants Ai, A2, and A3 determine the coupling strengths in the 
various spatial directions and allow to adjust the degree of anisotropy. It is often 
convenient to introduce a collective spin operator S with components 

1 ^ 
So.^i^Y.<^ «€ {1,2,3}, (4) 

i=l 

which allows us to rewrite the Hamiltonian ([T]) in the form 

Hh = -^ {\iSl + X2SI -f X^Sl) - 2hS^. (5) 

Certain choices for the coupling constants of the anisotropic Heisenberg model 
yield several important special cases, for example: (a) the isotropic Heisenberg model, 
Ai = A2 = A3; (b) the Ising model, Ai == = A2; (c) the isotropic Lipkin-Meshkov- 
Glick model, Ai = A2 and A3 = 0. For these special cases, and more generally whenever 
Ai = A2, the Hamiltonian ([5]) can be expressed entirely in terms of S'^ and S'3, i.e.. 



Nonequivalence of ensembles in the anisotropic quantum Heisenberg model 4 

the square and the S-component of the collective spin. As a consequence, an angular 
momentum eigenbasis diagonalizes H^ and S^ simultaneously, and the model can be 
solved by rather elementary means, as shown in section [3. II 

Here we consider the coupling constants Ai, A2, and A3 to be nonnegative, but 
otherwise arbitrary. The exact expression for the canonical Gibbs free energy g as a 
function of the inverse temperature /3 = 1/T|j and the magnetic field h is known for 
this model (and, in fact, for a larger class of systems) and is reported for example in 
[8] . The model is found to display a transition from a ferromagnetic to a paramagnetic 
phase in the canonical ensemble. 

3. Microcanonical entropy 

Owing to the long-range character of the interactions in the Hamiltonian ([1]), 
microcanonical and canonical ensembles cannot be expected to yield equivalent results. 
The main purpose of the present article is therefore to complement the known 
canonical results with those obtained in the microcanonical ensemble. 

In thermodynamics, the energy e is the variable conjugate to the inverse 
temperature /J, and the magnetization m is conjugate to —j3h. In the same way 
as (7(/?, h) represents the fundamental quantity of the quantum Heisenberg model in 
the canonical ensemble, the microcanonical entropy s(e, m) serves as a starting point 
for a microcanonical description in the thermodynamic limit. However, for a pair 
of variables (e,™) corresponding to the pair of noncommuting operators {Hq,M = 
Y] o'f), it is not even well established how to define a quantum microcanonical entropy. 
Extending a suggestion of Truong [5] to interacting systems, the definition 

SAr(e, m) = — In^ Tr [PHo{e)PM{m)] SA{e - e)SA{m - m) (6) 

e,m 

seems to be physically reasonable, but difficult to apply in practice. Here, e and 
m denote eigenvalues of the operators Hq/N and M/N, respectively. 6a is the 
characteristic function of the interval [— A, 0], i.e., Sa{x) = 1 if a; G [— A,0], and 
zero otherwise. The parameter A > is chosen small but otherwise arbitrary. Puaie), 
Phiipi) denote the eigenprojections of the operators Hq and M belonging to the 
eigenvalues e and to, respectively. 

Definition (j6|) can be motivated as follows: In classical physics, the density of 
states for a pair of variables is 

^classical (e, m) = '^ SA{e(x) - e)6A{m{x) - to), (7) 

X 

where the summation is over phase space (assumed here to be discrete for simplicity) 
and e and m are phase space functions. The interpretation of this classical expression 
is: Pick a state x with energy e{x) in the interval [e — A,e]. If furthermore 
fh{x) is in the interval [m — A,77i], then add one to the sum, otherwise zero. 
Examples of such densities of states (or the corresponding entropies Sciassicai(e,TO) = 
lnriciassicai(e, to)/A^) for classical spin systems can be found in [TUl lTniT^ . 

For a quantum system, the definition has to be modified: Since the operators Hq 
and M do not commute in general, an eigenstate of Hq is usually not an eigenstate 
of M. An eigenstate of Hq with eigenvalue e in the interval [e — A, e] will therefore 
in general not have a well-defined value of to. Instead, measuring the magnetization 

I Here and in the following Boltzmann's constant is set to unity. 
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of this energy eigenstate can yield one out of several values of to, each with a certain 
probability. It seems therefore reasonable to define the density of states [i.e. the 
argument of the logarithm in ([S])] by adding up, for all energy eigenstates with 
eigenvalue e in the interval [e — A,e], the probabilities for finding a magnetization 
in the interval [to — A, to]. For one-dimensional subspaces, this probability is given by 
the overlap |(e|TO)|^, and for higher dimensional subspaces (i.e. degenerate eigenvalues) 
it can be expressed in terms of projection operators as in the trace in (j6]). For the 
special case of commuting operators, these probabilities are either one or zero, and 
the density of states counts, in analogy to the classical case, the states that comply 
with a certain constraint. 

At least for finite systems, and in contrast to its classical counterpart, the entropy 
SM{e^ to) defined in ^ does not describe an ensemble of systems at fixed energy e and 
magnetization m: For noncommuting operators Hq/N and M/N, quantum mechanical 
uncertainty does not simultaneously allow fixed values for both quantities. It is, 
however, possible to give a probabilistic interpretation of the (suitably normalized) 
density of states, 

r2jv(e,TO) ex y^ Tr [Ph,, (e)PAf (to)] SaJb - e)SA{m- to), (8) 

as the probability of measuring, with experimental resolution A, a certain value of m 
at a given value of the energy e. 

In the thermodynamic limit N — > oo, the familiar interpretation of Sf^[e,m) as 
describing an ensemble of systems at fixed energy e and magnetization to, is recovered. 
One way to understand this property relies on an observation made by von Neumann 
in the early days of quantum mechanics ([13): see also |14j for a modern presentation 
of these ideas). Assume that an experimenter attempts to measure simultaneously the 
energy and the magnetization of a macroscopic system. Von Neumann argues that 
macroscopic observables do always commute, since nothing prevents the experimenter 
from reading out the display of his measurement devices simultaneously. However, the 
observables measured are not really energy and magnetization, but two commuting 
observables, approximating the observables of interest to a very good degree. De Roeck 
et al. |15j implemented this idea to define the microcanonical entropy of a quantum 
system in the thermodynamic limit, making use of the concept of concentrating 
sequences. In the same spirit, we would like to compute an entropy SAr(e, m), where 
e and to correspond to commuting operators Hq/N and M/N which converge to 
Hq/N and M/N in the thermodynamic limit in a suitable way. However, an explicit 
construction of Hq and M is difficult, and it is more convenient to use, as in ([6]), the 
observables Hq and M, and the resulting entropy coincides with the desired one in 
the thermodynamic limit. 

3.1. Special values of the coupling constants: Ai = A2 
By defining collective spin ladder operators 

5_ = i (5i - 152) , S+ = Sl = \{Si + iS2), (9) 

the Hamiltonian Hh in ([5]) can be written in the form 

Hh = -^ [(Al + A2) ^2 - (Al + A2 - 2A3) Si + 2 (Al - A2) [Si + Si)] ~ 2hS3. (10) 
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In the special case Ai = A2 = Aj^, the term containing S'_ and 5*+ vanishes, and the 
Hamiltonian can be expressed entirely in terms of the operators S'^ and S'3, 

Hh = -^ [X±S^ - {\i_ - A3) Si] - 2hS3. (11) 

The components Sa of the collective spin operator deserve their name, meaning that 
they obey the angular momentum algebra 

[Sa,Si3]=iea/3jSj, tt, /3, 7 G {1, 2, 3}, (12) 

as is easily verified from relation ^ and definition (|3]). Therefore, the well-known 
angular momentum eigenstates \S,M) are eigenstates of the Hamiltonian H^ with 
eigenvalues 

EsM,h = -|p [X±SiS + 1) - (Ai - A3) m2] - 2hM. (13) 

From the rules for the addition of angular momenta it follows that S can take on 
integer values from to N/2 where, for simplicity, we have assumed N to be even. 
The Hilbert space can be decomposed into a direct sum of irreducible representations 
of S'C'(3), 

N/2 

H = (C2)®^-0dsI?s. (14) 

s=o 
Each representation Vg is a linear subspace of dimension 25" + 1 with basis \S,AI), 
Me {-S.,...,+S}, and 

m{2S+l) 

^ {N/2 + s + iy.{N/2-syr ^ ' 

is the multiplicity of the representation Vs for N even [T^ . 

Now the trace in ^ can be written with respect to the \S^ M )-eigenbasis, 

N/2 +s 

f]jv(e,™) = ^ds Yl ^A{Ne~EsM,Q)SA{Nm~2M). (16) 

In the limit of small A and large system size N, we can write 
s{e,m) = lim — ln51jv(e, to) 

N^oo N 

- hm 1 m y ds S. (Ne + ''-'^' + '^ ^"^^ ^'^ ' '^^ 



1 , N12S 

hm — In ■ 



w^oo N {N/2 + S)l{N/2 - Sy. 



2S=Ny/m'^{l-X3/Xj_)-2e/\^ 

Applying Stirling's formula, we obtain for the microcanonical entropy of the aniso- 
tropic quantum Heisenberg model in the thermodynamic limit the final expression 

s(e, to) = ln2 - i[l - /(e, m)] ln[l - /(e, to)] - ^[1 + /(e, m)] ln[l + /(e, to)] (18) 
with 



where Aj^ = Ai = A2. The domain of s is given by 

P = {(e, m) e R^ I 2e -I- m^A3 < and A^ > to^ (A± - A3) - 2e} . (20) 



Nonequivalence of ensembles in the anisotropic quantum Heisenberg model 7 

3.2. General values of the coupling constants 

For positive, but otherwise arbitrary, values of the couphng constants A^ in the 
Hamiltonian ([T]), an evaluation of the expression © is difficult to achieve. Considering 
the thermodynamic limit of the entropy, 

s{e,m) — lim sjv(e,?Ti), (2f) 

usually accounts for a simplification of the problem, but not enough so in this case. 
Here, in order to render the calculation feasible, we claim that the exact result of 
expression pi|) is recovered by formally computing, in a sense which will become clear 
in the following, the quantity 

s(e,m) — lim — In J7jv(e, to) (22) 

with a density of states J7jv that we write symbolically as 

rjjv(e, to) = Tr [5{e - Ho/N)d{m - M/N)] . (23) 

Note that the symbolic expressions make little mathematical sense and require some 
physically reasonable regularization, like the A-regularization in ([8]) . Here, instead, we 
will keep the unregularized deltas for the moment to render some formal manipulations 
more obvious. Expressing the deltas in their Fourier representation and performing 
the thermodynamic limit later on will cure the problem and lead to a well-defined 
mathematical expression for the microcanonical entropy s(e, to,) in the thermodynamic 
limit. 

In view of the deltas, the integral remains unchanged upon insertion of the factors 
exp[a(A^e — Hq)] and exp[&(7VTO — M)] in the trace, 

flNie, to) = Tr {N^ exp[a(iVe - Ho)]S{Ne - Hq) exp[b{Nm - M)]S{Nm - A/)} , (24) 

where a and b are real parameters. Writing the deltas in Fourier representation, we 
obtain 

r2 



ilNie, m) = ^ f dk f di Trjexp (a + ifc) ( A^e + ^ {Sf + 5| + 5|) j 



X exp [{b + i£) (Nm - 2S3)] 



(25) 



where the anisotropic collective spin operators 

Si = A,^^ (26) 

have been introduced to ease the notation. Unless specified otherwise, domains of 
integration always extend from —00 to +00. 

3.2.1. Decoupling the N-spin trace. To evaluate this expression, we want to rewrite 
the density of states ()25|) in such a way that the A^-spin trace decouples into a product 
of one-spin traces. To this purpose, we can adapt to our microcanonical setting a 
couple of tricks that have been used by Tindemans and Capel [171 in a related canonical 
calculation. The first exponential in (|25p is an exponential of a sum of squares of the 
collective spin components Sa, which — for a reason that will become clear soon — we 
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would like to rewrite as a product of exponential operators. This can be done by using 
a straightforward generalization of the Lie- Trotter formula [IH] , yielding 



exp 



{a + \k)[Ne+-[St+S'^+Sl 



— lira 



exp 



(a + ifc)iVe 



n 

Q = l 



exp 



2(a + ifc) 
nN 



St. 



(27) 



Now the exponential operators are in a suitable form for applying the Hubbard- 
Stratonovich trick [191 HO] ■ This is a fancy name for the identity 

1 



e^o = 



dxe-" /(™ )e2-o/", 



(28) 



where O is an operator and c is a constant with 5R(c) > 0. With this formula, we can 
express ((27)) as 



lim 



N 



2Tm{a + ifc) 



3/2 



^{a+ik)Ne/n 



a; exp 



Nx ■ X 
2n{a + ik) 



n 



,2x^Sc/r, 



, (29) 



where x — {xi,X2,X3). As a result, the operators Sa now occur linearly in the 
exponent, which will be important for what follows. Next we replace the product 
(over a) of three exponentials in (1^^ by the single exponential 



exp - ^ XaSa 



(30) 



This term is in general not identical to the product in ()29|) . but the correction — 
as given by the Baker-Campbell-HausdorfF formula — is of the order 1/n^ and should 
therefore in the limit n — >■ oo be negligible when compared to terms of the order 1/n. 
A rigorous justification of the replacement (|30)) is given in Appendix A of [17] . 
Recalling that the operators 2Sa = 



^N 



operators a", we can rewrite the exponential pop as 

exp 
where 

3 
^{i,x) = -^Xa\/)^(J°' 



\ a=l i=l / \ i=l 



Aq X]i=i ^f ^'^^ sums of one-particle- 
e(j,a;)), (31) 



(32) 



Since the ^(i, x) are one-spin operators, acting non-trivially on only one of the factors 
of the tensor-product Hilbert space ((D^)*^^, their commutators must vanish, 

K(«,a;),C(j,y)] =0 io^i^j. 

Hence, the exponential in pip is equal to 

N 



n< 



and we obtain 



-i(i,x)/n 



N 



(33) 
(34) 



27rn(a -I- ifc) 



3/2 



d xexp 



Nx ■ X 
2n{a + ifc) 



N 



-^{i.x)/n 



(35) 
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for the expression in (1291) . Now we can express the nth power of a three-dimensional 
integral in (I55|) as a Sn-dimensional integral, 

3n/2 

(36) 



n^oo y 27rn{a + ik) 



(a+ifc) ;ii' 



2n(c 

Inserting this expression into the density of states psp and making use of the linearity 
of the trace, we can write 



fiAr(e,m) = — . / dse^'^'^ / die^™* lim < 

47r^ Ja-ioo Jb-ioo "-^°° 



N 



2'Kns 



3n/2 



• • • / dVi) • • • dV") 



iV 



exp 



One -"^ 



xTr 



Cn N \ 

m=li=l / 



2ns 



2tS3 



,(™) . a.(™) 



(37) 



Here we have also substituted the integration variables k and ^ by s = a + ifc and 
t = b + i£. The trace in ((37)) acts only on exponentials of one-particle operators, and 
this allows us to rewrite the trace Tr on the iV-spin Hilbert space ?^ as a product of 
traces tr^ over one-spin Hilbert spaces, 



Tr 



N 



N 



]J]-Je C(,x<">)/« Y[e 



\7n—l i—1 



N 



Tri J] ne-[«(^'-""')+*-' 



lA 



\ni—l i—1 



N / n \ 

UtvA J] e-[«(''^""')+*-.']/" 



(38) 



N ( r n / 3 

ntr. -p ^E E-^'^ 

Z— 1 l^ L 771 — 1 \a — 1 



' V Aa (7, — 1(7, 



For the same kind of reasoning as outlined below pOI) . commutators of order 1/n^ 
have been neglected. In the limit n — >■ oo this approximation will become exact and 
the asymptotic equalities ~ in psp become proper equalities. With the definition 



(39) 



1 , " 



we can rewrite (|38|) as 

N p 

exp 



ntr 

z=l 

where 



C3 Ci - 1C2 

Cl + iC2 -C3 



tr 



N 



N 



(2coshrt) , 



n^n{{x("^^})^ ^4 + 4 + 4 



(40) 



(41) 
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Replacing the trace in ([57]) by expression (1^0]) , we obtain for the density of states 

(42) 



oAf 1^2 /■a+'ioc pb+ioo 

Vli^{e,m) = ^ / As j At 

J a — ioo Jb—ioo 



X lim f ^^] " / • • • / d^x^i) • • • d^x^n)^N^isM<^ 



where 

n 

JF{s,t,{x("''>}) = es + mt-—- V x^""^ ■ x^""^ +\ncosh\rt(\x'^"''>])]. (43) 

m— 1 

3.2.2. Asymptotic evaluation of the (3n + 2) -dimensional integral. We have to deal 
with a (3n + 2)-diniensional integral, to be evaluated in the limit n — > oo and, 
afterwards, the thermodynamic limit N —> cx). This integral (|42]) is of Laplace- 
type with respect to the large parameter iV, and an asymptotic evaluation can be 
performed by a multidimensional version of the method of steepest descent (see |21] 
for a textbook presentation). 

To apply this method, we need to find a stationary point of the function T 
for which it is possible to smoothly deform the contours of the s- and i-integrations 
such that the paths of integration correspond to constant (zero) imaginary part of J-. 
Stationary points of J- need to satisfy the conditions 



= —— = e 



ds 2ns 



n 

— Vaj^^^-a;^"), (44a) 



m—l 



dF tanhrf / y/X^ ^ (,„) 

at rt \ n 



E 4™^ - (44&) 



m— 1 



dF Xa tanhrt "v/A^ / ^fKi 



dXa ^^ Tt n \ n ^-^ I 



771 — 1 



where a G {1,2,3} and u E {1, . . . ,n}. Inserting (144 &P in (|44c|) and some straightfor- 
ward algebra allows us to rewrite this set of equations as 

n 

= 2nes^ + ^ a;'™' • a;^"'\ (45a) 



771—1 



a 



= mri+ t-^-^E4"Mtanhrt, (456) 



n 

771 — 1 



= x^"' - msy%, (45 c) 

= 4"' (^ - ^ E 4"^1 + ^ E 4-^ (45rf) 



n 

771 — 1 



where a G {1,2} and m G {1, . . . , n}. There is a class of particularly simple solutions 
to this set of equations where all the a;*^™' are identical, i.e. 

x^"^^ ^x^{xi,X2,X3) VmG {l,...,n}. (46) 

Similar to the canonical calculation reported in jT7], it should be possible to prove 
that one of the stationary points of T which is subject to condition (|^5)) corresponds 
indeed to the maximum of the exponent along the properly deformed integration path 
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and therefore yields the correct resuh for the integral (j42p in the thermodynamic limit. 
We have not gone through this calculation explicitely, but the final result, especially 
when compared to the results of section [331 seems to confirm this assumption beyond 
any reasonable doubt. Under the assumption ([^S)). the set of equations simplifies to 

= 2es2 + a;^ (47a) 

= mRt{x) + (t- y^xs] ta.nhRt{x), {47b) 



where 



= 2^3 — msy'X:}, (47c) 

= Xa[ms{X3~Xa)~t], ae{l,2}, (47rf) 



Rt{x) ^ ^ Xixj + X2xl + (^t - ^3x3) . (48) 

For an asymptotic evaluation of the integrals in (|42p by means of the method of 
steepest descent, we have to evaluate J"(s,i, {aj^"*)}) as defined in (|43|) at the values 
of s, t, and {x^™^} specified by (|47ap - (|47(i|) - For J^ we obtain under condition (|46| 
the expression 

J^{s, t,x) = es + mt- In [1 - tanh^ Rt (x)] , (49) 

where the identity 2 In cosh x — — ln(l — tanh x) has been used. Making use of the 
identities (|47a|) - (j47rfj) . it is a matter of straightforward algebra to evaluate T at the 
values So, to, and Xq which are solutions of this set of equations (see [Appendix A| for 
the details of the calculation) . Evaluating J^ at these points and taking into account 
(|T7at - ([T7c[) . we obtain 

- 2J-(so, to, xo) = [1 - fie, m)] ln[l - f{e, m)] + [1 + /(e, m)] ln[l + /(e, m)] (50) 

with 



/(e,m) = Jm2(l-^)-^ (51) 



as defined previously in (1191) . The constant Aj^ can be either Ai or A2, unless one 
of them vanishes. The solution relevant for the asymptotic evaluation of the integral 
(|42p is the one which maximizes J^. It is straightforward to verify that this solution 
corresponds tq§| 

Ax = maxjAi, A2}. (52) 

Solutions of the type ([50l) exist for all values of {e,m) which satisfy the inequalities 

2e + m^As < and A± > m^ (Ax - A3) - 2e (53) 

(see again [Appendix~A[ for a derivation). 

According to the method of steepest descent, the asymptotic behaviour of 51 jv in 
(|42p is now given as exp[A^(ln2 + J^)] times some prefactor (see for example section 
3.7 of Miller's textbook [21] for the prefactor of multidimensional Laplace integrals, 
which can be adapted to the method of steepest descent of multidimensional integrals). 

§ The case of Ai = = A2 has to be treated separately, but this is just the well-known case of the 
Curie- Weiss Ising model. 
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The prefactor, however, is subexponential in N . Since we are interested in the 
microcanonical entropy in the thermodynamic hmit, 

s{e,m) = hm — - Inrijvfe, m), (54) 

subexponential terms do not contribute, and we obtain the following final result 
for the microcanonical entropy of the anisotropic quantum Heisenberg model in the 
thermodynamic limit; 

s(e,m)=ln2-i[l-/(e,m)]ln[l-/(e,TO)]-i[l + /(e,m)]ln[l + /(e,m)], (55) 

defined on the domain 

P = { (e, m) G r2 I 2e + m^ A3 < and A^ > m^ (A_l - A3) - 2e} , (56) 

where /(e, m) and A_l are specified in ([CTI) and (15^ . 

This result is remarkably simple, in the sense that an explicit expression for s(e, m) 
can be given. This is in contrast to the canonical ensemble, where the canonical free 
energy g{/3, h) is given implicitly as the solution of a maximization (see section 2] 
for more details on the canonical solution). It may seem a bit disappointing that 
the solution (j55p is identical to the one we obtained already from the much simpler 
calculation for the special case Ai — A2 in section [O] with the only difference that 
now Aj_ is defined according to equation ([52]) . Plots of the domains and graphs of 
s{e,m) are shown in figure [I]for a number of coupling strengths A^, A3. 

4. Recovering the canonical Gibbs free energy 

The microcanonical entropy we have computed in section |3] forms the starting point 
of an analysis of the anisotropic quantum Heisenberg model in the microcanonical 
ensemble. An analogous role is played by the canonical Gibbs free energy 

gN{l3,h) = -^InTre-^"- (57) 

for calculations in the canonical ensemble. In the thermodynamic limit, the 
corresponding infinite-system quantity 

5(/3,/i)= hm gNiP,h) (58) 

A/— ^00 

is related to the microcanonical entropy s{e,m) by means of a Legendre-Fenchel 
transform, 

- Pg{P, h) = sup [s(e, m) - /3e + /3/ito] . (59) 

e.m 

Note that while, in the thermodynamic limit, the canonical free energy is always given 
as the Legendre transform of the microcanonical entropy, the inverse is in general not 
true. In particular, it would have been impossible to derive the entropy ([55)) from the 
canonical Gibbs free energy for anisotropy parameters A_l < A3 for which the entropy 
is nonconcave. This will be discussed in more detail in section [5] 

In order to determine g, it is helpful to consider the following operational 
interpretation of the supremum in (|59p : Take a plane described by the equation 

Si3,h,c{e,m) — c + l3e — I3hm, (60) 

i.e. with slopes /3 in e-direction and —/3h in m-direction. Start with a very large value 
of the parameter c and lower this value, and therefore the plane, until, at some value 



Nonequivalence of ensembles in the anisotropic quantum Heisenberg model 



13 




-0.5-0.4-0.3-0.2-0.1 0.0 




-0.5-0.4-0.3-0.2-0.1 0.0 




-0.5-0.4-0.3-0.2-0.1 0.0 




-0.5-0.4-0.3-0.2-0.1 0.0 




-0.5-0.4-0.3-0.2-0.1 0.0 




-0.5-0.4-0.3-0.2-0.1 0.0 










-0.5-0.4-0.3-0.2-0.1 0.0 




Figure 1. Domains X? (left) and graphs (right) of the microcanonical entropy 
s{e, m) of the quantum Heisenberg model for some combinations of the couphng 
Ax, A3. From top to bottom; (Ax, A3) = (1/4,1), (9/10,1), (1,1), (1,9/10), 
(1,1/2), (1,1/5), (1,0). For the domains, the abscissa is the energy e and the 
ordinate is the magnetization m, and the entropy is defined on the shaded area. 
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Figure 2. Plot of the boundary curve described by \62) and l|63ll together with 
two examples of tangent planes touching the curve at the point (e(mo), "JOi *("io)) 
with mo = 1/2- The planes have slopes /3 = 1/5, respectively /3 = 1/2, in the 
e-direction, while their slopes in m-direction are determined by II66I I. 



of c(/3, /i), the plane touches the graph of s(e, m) for the first time. Then the value of 
g{P, h) is related to Sjjhc{i3,h) evaluated at the origin, 

-/3g(/3,/l)=S^,;,,e(/3,/^)(0,0). (61) 

More details on this graphic-geometric interpretation, although for the special case 
of a Legendre (not Legendre-Fenchel) transform, can be found in section III. A and 
figure 3 of [22] . We will now use this interpretation to calculate the Legendre-Fenchel 
transform of s, separately for the cases A_l < A3 and Ax > A3. 



4.1. \i_ < A3 

By inspection of rows one and two of figure [T] [or by analysis of the results in ()55|) 
and ([56]) ] ■ one can convince oneself that, for X± < A3, a plane of the type (|60l) when 
lowered onto the graph of s will always touch the graph at a boundary point of the 
domain I? of s, given by 

e{m) = -^m^Xs, (62) 

where 

s{m) = s (— iTO^A3,m) = ln2 — i(l — m) ln(l — m) — i(l + m) ln(l -|- m). (63) 

Next we need to construct the family of all planes tangent to this boundary curve at a 
given point m = toq; see figure [5] for an illustration. A tangent vector to the boundary 
curve is 

t{mo) — edmoe{mQ) + rhdmotTiQ + J9,„Qs(mo) = — empAs +rh— sarctanh(r7T,o), (64) 

where e, m, and s are unit vectors in e, m, and s direction. A family of vectors 'n{mo) 
which are perpendicular to t{mQ) can be constructed by defining 

n{mo) = —f3e + (3hm + s (65) 

and requiring that 

n(mo) • t(mo) = /3moA3 + (ih — arctanh(TOo) = 0. (66) 

The constraint ([55)) fixes one of the parameters /3 and h in ([55)) . while the other 
parameter labels a family of (non-normalized) vectors. To obtain the family of tangent 
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Figure 3. Canonical phase diagram of the Curie-Weiss anisotropic quantum 
Heisenberg model. The lines in the (rescaled) {/3, h)-planes indicate the values of 
inverse temperature and magnetic field at which the canonical Gibbs free energy 
g(/9, /i) is nonanalytic. Left: Phase diagram for anisotropy parameters Aj^ < A3. 
For small values of /3A3 , the system is in a paramagnetic phase, whereas it becomes 
ferromagnetic for larger values. Right: Phase diagram for anisotropy parameters 
-^± > ^3- Inside the C-shaped curve, a non-zero magnetization in the 1-2-plane 
occurs, whereas outside this area the magnetization vector points into the 3- 
direction. 



planes to the boundary curve, we consider each of the vectors n{mo) as a normal vector 
of such a plane. We further demand that each plane touches the curve in the point 
p(mo) = (e(mo), mo, s(mo)). All points r = (e, m, s) lying in this plane have to fulfill 
the equation 

= n(mo) • (r - p(too)) = -(3 {e+ ^mlXs) + (3h{m - mg) + s - s{mo). (67) 

Therefore, the planes s touching the graph of s as explained above are given by 

Si3,h = s(mo) + /3 (e + ^mlXs) - (3h{m - mo), (68) 

where /3, h, and m,o are subject to (|66|) . For /3A3 < 1 or for sufficiently large magnetic 
fields h, equation ([66|) has just a single solution, whereas for PX^ > 1 and small enough 
h three solutions exist. In the latter case, the relevant plane s corresponding to the 
supremum in (j59l) is given by the mo with the largest absolute value. According to 
(|6ip we have 

- /3g(/3, h) ^ s(mo) + i/^m^Ag + l3hmo (69) 
and, making use of (j66p . this result can be written as 

- /3g(/3, /i) = ln2 - i In (1 - m^) - i/Jm^Ag, (70) 

where mo = rnQ{/3, h) is determined implicitly as the one solution of (j66p having the 
largest absolute value. The swapping from one solution branch to another at ft, = 
for inverse temperatures /3A3 > 1 gives rise to nonanalytic behaviour of the canonical 
Gibbs free energy density, resulting in the phase diagram as plotted in figure [3] (left). 

Note that, in contrast to the explicit expression for the microcanonical entropy 
([551) . the transcendental equation ([M)) cannot be solved explicitely for mo, and 
therefore only an implicit expression for g{j3, h) can be given. 



4-2. Ai > A3 

By inspection of rows four to seven in figure [U one can infer that, for anisotropy 
parameters A_l > A3, there is a certain range of slopes — /3 and /3h for which a plane 
of type (|60p touches the graph of s(e, m) not at a boundary point of the domain V of 
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s as in section [4.11 but in its interior. This is true for precisely those values of j3 and 
h for which the equations 

ds{e,m) arctanh[/(e,TO)] 
de Xj_f{e,m) 

ds{e,m) /As \ arctanh[/(e, m)] 

dm \X± J f[e,m) 

have sokitions in V. This set of equations can be rewritten as 

h={\j_~ A3)m, (72a) 




- 2e ) = arctanhJ -— ( ^ 2e ) , (72&) 

A3 / y Aj^ \A^ - A3 / 

and one can verify that, for Aj_ > A3, the equations have solutions with {e,m) G V 
for all /3 > and h satisfying 

P ^ ^T^ arctanh ( -^] ■ (73) 

hXi_ \A_L-A3/ 

For these values of /3 and h, the Legendre-Fenchel transform ((59|) reduces to the 
more familiar Legendre transform 

- /?5(/3, h) = s(e(/3, h), m(/3, h)) - /3e(/?, h) + fihm{(i, h), (74) 

where e{(3,h) and m{j3,h) are solutions of equations (|72ap and (|72&p . Rewriting the 
entropy (j55p in the form 

s(e, to) = In 2 + /(e, to) arctanh [/(e, to)] - i In [l - /(e, m)^] (75) 

and making use of (|72ap and p2b\ , we obtain as a final result for the canonical Gibbs 
free energy 



/3g(/3,/i)=/3e(/3,/i)~iln 



A / h 



2 



4-aiIa7^-2^(^'^) 



(76) 



where e{(3,h) is a solution of (J72&|1 . Equation ([76l) is valid for all (/?, /i) satisfying 
the inequality (|73p . For other values of inverse temperature and magnetic field, it is 
again the boundary ([62|) of the graph of s(e, to) that determines the Legendre-Fenchel 
transform, and the result for g{/3, h) is the one derived in section 14.11 and stated in 
(iTOll. The curve 



/3/iA_L = iX± - A3) arctanh — (77) 

V Aj. - A3 / 

in the (/3, /i)-plane, at which the inequality (j73p becomes sharp, determines the 
boundary separating the two different kinds of behaviour of g. This is precisely the line 
in the phase diagram separating the ferromagnetic phase of the Heisenberg magnet 
from the paramagnetic one, as plotted in figure |3] (right). 

5. Nonequivalence of ensembles 

In the previous section, we derived the canonical Gibbs free energy g(/3, h) from the 
microcanonical entropy s{e,m). While, in the thermodynamic limit, it is always 
possible to compute canonical thermodynamic potentials from microcanonical ones, 
the inverse is not necessarily true. In those cases where the backward transition (from 
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canonical to microcanonical) is also possible, one speaks of equivalence of ensembles, 
otherwise of nonequivalence. 

By inspection of rows three to seven in figure [T] [or by simple analysis of the 
results in ([55t and ([55]) ] . the entropy s for A^ > A3 is seen to be a concave function 
on a domain which is a convex set. For A^ < A3, the domain is not a convex set and 
therefore the entropy is neither convex nor concave. In the latter case, microcanonical 
and canonical ensembles are not equivalent, in the sense that it is impossible to obtain 
the microcanonical entropy s(e,m) from the canonical Gibbs free energy g{l3,h) by 
means of a Legendre-Fenchel transform. This is evident from the fact that the outcome 
of a Legendre-Fenchel transform is always concave or convex [23] . 

The physical interpretation of ensemble equivalence is that every thermodynamic 
equilibrium state of the system that can be probed by fixing certain values of e and 
m can also be probed by fixing the corresponding values of the inverse temperature 
(3{e, m) and the magnetic field h{e, m). In the situation A_l < A3, where nonequivalence 
holds, this is not the case: only equilibrium states corresponding to values of (e, m) for 
which s coincides with its concave envelope can be probed by fixing (/3, /i); macrostates 
corresponding to other values of (e, m), however, are not accessible as thermodynamic 
equilibrium states when controlling temperature and magnetic field in the canonical 
ensemble. In this sense, microcanonical thermodynamics can be considered not only 
as different from its canonical counterpart, but also as richer, allowing to probe 
equilibrium states of matter which are otherwise inaccessible. For more information 
on nonequivalent ensembles, see [24j or the introductory article [25j . 

There is a further long-range peculiarity, going under the name of partial 
equivalence [26j , which can be observed in the anisotropic quantum Heisenberg model 
for any values of the coupling constants Ai, A2, and A3. Partial equivalence here refers 
to the situation where a macrostate, associated with a certain pair of values (e, m) 
in the microcanonical ensemble, corresponds to more than just one pair of values 
(^, h) canonically. We have encountered this situation repeatedly when calculating 
the canonical Gibbs free energy g. In particular for the case Aj^ < A3 in section WA\ 
we computed "touching planes" s with different slopes /3 and —fih, all touching the 
graph of the entropy s at the same boundary point (— TOgA3/2, toq). As a consequence, 
all pairs (/3, h) which solve equation ([66]) for the same value of mo correspond to the 
same microcanonical macrostate labelled by the parameters (e(mo)j"T'o)- The same 
holds true for A^ > A3 and values of (3 and h for which 

^,^„U„l,(^). ,78) 

This provides also an explanation for a peculiar observation reported in a study of 
the fidelity-metric of the Lipkin-Meshkov-Glick model, i.e. the special case A3 = 
[27] . This article reports that, for /? and h satisfying (|78| . a fidelity metric on the 
(/3, /i)-plane is not well defined. Having observed partial equivalence in this parameter 
regime, we know that entire curves, not single points, in the (/3, ft,)-plane label one and 
the same macrostate. As a consequence, all (identical) macrostates along such a curve 
have distance zero from each other, which implies that a metric on the (/3, /i)-plane 
cannot be positive definite and is therefore not well defined. 
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6. Thermodynamic equivalence of Heisenberg and Ising models 

It had been observed already in the 1970s that the Curie- Weiss isotropic Heisenberg 
model (Ai — A2 = A3) and the Curie- Weiss Ising model (Ai — = A2) are 
thermodynamically equivalent in the sense that their canonical Gibbs free energies 
coincide in the thermodynamic limit [28]. From the results we have obtained in 
section |4l it is obvious that an even more general statement can be made: For all 
anisotropy parameters satisfying Aj^ < A3, the canonical Gibbs free energy g{f3,h) 
of the Curie- Weiss anisotropic Heisenberg model coincides with that of the Ising 
model. This is a direct consequence of the fact that, as pointed out in section 14. 1[ 
the Legendre-Fenchel transform of s{e,m) for Aj_ < A3 is determined exclusively by 
s(— 171^X3/2,171), i.e. the entropy evaluated at the boundary e(m) = —Tn^XzjI. Since 
the entropy at this boundary is identical to the entropy of the Curie- Weiss Ising 
model, thermodynamic equivalence holds in the canonical ensemble. Putting this 
differently, the microcanonical entropies s(e, to) of the Curie- Weiss Ising model and the 
Curie- Weiss Heisenberg model with A^ < A3 share identical concave envelopes, which 
implies that their Legendre-Fenchel transforms must be the same [23] • Remarkably, 
however, thermodynamic equivalence does not hold in the microcanonical ensemble, 
as is obvious from the different shapes of entropies in rows one to three of figure [TJ 

Vertogen and de Vries have claimed in reference '291 that the Curie- Weiss Ising 
model and the anisotropic Curie- Weiss quantum Heisenberg model are equivalent not 
only thermodynamically, but even on the level of their Hamiltonian operators in the 
thermodynamic limit. It is not obvious to the author how to reconcile this claim 
with the results of the present article: Equivalence on the operator level should imply 
identical thermodynamic functions also in the microcanonical ensemble, in conflict 
with the results of section |3l 

7. Physical relevance and experimental realization 

In the thermodynamic limit, the microcanonical entropies ([5]) and (j22p discussed in 
this article describe the physical situation of fixed energy e and fixed magnetization to,. 
For neither of the two constraints it is immediately obvious how they can be realized 
in experiment: The quantum Heisenberg model was devised to model ferromagnetic 
spin systems which, in their traditional condensed matter realizations, are typically 
coupled to a thermal reservoir. As a consequence, the energy is not fixed, but fiuctuates 
around a certain mean value, and the canonical ensemble is appropriate for a statistical 
equilibrium description of this situation. 

Recently, however, it has been pointed out that cold atoms in optical lattices 
are an ideal laboratory for engineering systems which are governed by Hamiltonian 
operators formally equivalent to those of condensed matter spin systems. Furthermore, 
such cold atom realizations of condensed matter-type systems possess the very 
attractive feature of being highly controllable: by appropriately tuning Feshbach 
resonances and other parameters, the interaction type and strength can be tuned freely, 
even changing the character of the interaction force from attractive to repulsive. After 
switching off the cooling in such an experiment, total energy and number of atoms are 
conserved to a very good degree. As a consequence, a statistical description of such a 
lattice spin model should make use of the microcanonical ensemble. 

The forces between cold atoms in optical lattices are typically of very short range: 
In most cases, s-wave scattering is dominant, and in this case it is appropriate to 
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model the interaction force by an effective contact interaction. As a consequence, 
also the total magnetization is a conserved quantity, and the ensemble of constant 
energy and magnetization is appropriate for describing the equilibrium properties of 
such systems. A discussion of short-range spin systems under these conditions can be 
found in [301 [12]. 

Nonequivalence of ensembles is however restricted, as explained in the Introduc- 
tion, to long-range interacting systems. A long-range interaction can be introduced 
in cold gases by using atoms or molecules with a permanent electric or magnetic 
dipole moment, resulting in a dipole-dipole interaction potential decaying like r~^ 
with the interparticle distance r. Alternatively, as shown by O'Dell et al. |31| . an 
r~^-interaction can be engineered by inducing, with appropriately tuned laser light, a 
dipole moment in atoms without a permanent dipole moment. Although such an 
algebraic decay is obviously different from the distance-independent forces in the 
Curie- Weiss-type model we have studied, it is known that Curie- Weiss-type models 
faithfully reproduce many properties of algebraically decaying long-range interactions 
qualitatively, and to some extent even quantitatively |321 133] . In particular, one 
can show under rather mild conditions that, by making the parameter a > of 
an algebraically decaying r~"-interaction small. Curie- Weiss behaviour is approached 
continuously. Therefore, if the Curie- Weiss model shows a nonconcave entropy, the 
same is to be expected for algebraically decaying interactions with sufficiently small 
values of a. Of course, this does not guarantee that nonequivalence of ensembles 
indeed persists to a = 1 or even a = 3, but it appears to be at least a plausible 
scenario. 

It has been shown by Micheli et al. ^ that dipolar gases in optical traps allow for 
the realization of many Hamiltonians of interest in condensed matter physics, including 
the anisotropic quantum Heisenberg model. Experimentally, impressive progress in 
cooling and trapping cold dipolar atoms and molecules has recently been made |34j |35] 
and, although dipolar gases in optical lattices have not been realized at the time 
of writing, it seems just a matter of time until this goal will be achieved and the 
engineering of long-range interacting anisotropic quantum spin models will be possible. 
Unfortunately for our purposes, the dipole-dipole interactions allow for scattering 
transitions beyond s-wave scattering, and the total magnetization is therefore not a 
conserved quantity. As a consequence, the statistical ensemble realized in such an 
experimental setting is a microcanonical one with constant energy, but fluctuating 
magnetization. The entropy 

s(e) = maxs(e, m) (79) 

ra 

appropriate for describing this situation is, however, a concave function, and 
nonequivalence of ensembles cannot be detected under such conditions. 

To probe nonequivalence in a cold atom experiment, we need to resort to a long- 
range spin model showing a temperature-driven discontinuous phase transition: In this 
case, the microcanonical entropy s(e) can be expected to be a nonconcave function, 
and an observation of this property would be possible under the experimental condition 
of conserved energy, but fluctuating magnetization. Such models will be investigated 
in a future work. 
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8. Conclusions 

The main goal of tlie present article was to contribute towards the understanding 
of nonequivalence of ensembles in quantum spin systems. To this purpose, an 
exact, analytic calculation of the microcanonical entropy s{e,m) of the anisotropic 
Curie- Weiss quantum Heisenberg model in the thermodynamic limit was reported. 
To the best of the author's knowledge, this is the first microcanonical calculation 
reported for a quantum spin system. With slight modifications, the strategy used 
to solve this problem, partly inspired by a canonical calculation by Tindemans and 
Capel [T7], should also be applicable to other quantum systems with Curie- Weiss- 
type interactions, for example to Heisenberg models with spins larger than 1/2. 
From the microcanonical result, the canonical Gibbs free energy (;(/3, h) was then 
recovered by means of a Legendre-Fenchel transform. This transform, especially in 
those instances when it does not coincide with a Legendre transform, also provides 
the key to understanding certain kinds of equivalences or nonequivalences that show 
up in the anisotropic Curie- Weiss quantum Heisenberg model: We found that the 
microcanonical entropy s{e,m) is a nonconcave function for anisotropy parameters 
A± < A3, and in this case microcanonical and canonical ensembles are nonequivalent. 
Furthermore, independently of the values of the anisotropy parameters, partial 
equivalence occurs: different pairs (/3, /i) and (/3',/i') of the canonical variables 
can correspond to the same pair (e, to) of microcanonical variables. Finally, the 
microcanonical result sheds light on the observation, dating from the 1960s and 1970s, 
that the Curie- Weiss Ising and the Curie- Weiss quantum Heisenberg models have 
identical canonical Gibbs free energies. We have shown that this thermodynamic 
equivalence also holds true in the anisotropic case for anisotropy parameters A^ < 
A3, and it is a consequence of the peculiar nonconcave shape of s(e, tti) for these 
parameter values. Microcanonically, however, Heisenberg and Ising models are 
not thermodynamically equivalent, as is obvious from their differing microcanonical 
solution ([55]) . 

Microcanonical solutions of quantum spin model are argued to be relevant for 
a statistical description of dipolar gases in optical lattices. In such experiments, 
the energy is controlled and conserved to a very high degree, rendering a micro- 
canonical description adequate. However, the peculiarities of long-range systems, like 
nonconcave entropies, nonequivalence of statistical ensembles, or negative microcan- 
onical response functions, should not depend on the precise nature of the long-range 
interactions. Still, a study of nonequivalent ensembles in quantum spin systems with 
algebraically decaying long-range interactions, as potentially realized in optical lattice 
experiments, is of course worthwhile and planned for future work. 

These results and discussions point out the importance of nonstandard thermo- 
dynamics beyond the canonical ensemble for experiments with cold dipolar atoms 
or molecules in optical lattices: Equivalence of ensembles does not hold in general, 
and a comparison of experimental data with canonical statistical physical predictions 
is bound to fail in this case. On the other hand, the results show that such cold 
atom experiments can provide an ideal laboratory for studying fundamental issues of 
thermostatistics, like nonequivalence of ensembles, in a highly controlled environment. 
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Appendix A. Evaluation of T at the stationary points 

It is shown how to evaluate J- as given in (|49|) at a stationary point determined by 
equations (|47ap - (|47(i|) - There are several ways to satisfy ()47(ij) for a = 1, 2; 

(i) xi — Q and ms{\^ — A2) = i, 

(ii) X2 — Q and nis{\^ — Ai) = i, 

(iii) 7Tis(A3 — Ai) — t and tos(A3 — A2) = t, 

(iv) xi — and X2 — 0. 

We discuss (i) in detail and argue later that the other cases do not contribute anything 
new. 

Assuming that xi — and tos(A3 — A2) = t, and making use of ()47c|) . equation 
(|ig|) simplifies to 



Rt = ^\2{xl+m-^ 3-^X2), (A.l) 

and the set of equations (J47al) - (j47(ij) takes on the form 

= 2es2 + xl + xl = s'^ (2e + m^As) + xl, (A.2) 

0==m(i?t-sA2tanhi?t). (A.3) 

From (|A.2p it follows that solutions exist only under the condition 

2e + m^Aa < 0. (A.4) 

Inserting (|A.2I) into (lA.ll) we obtain 



^^H-^-l^ 



s\2 



(A.6) 



which ahows us to write J^, given in (|49|) . in the form 

T = es + mt ^--In ' ' 

2s 2 

— —es 

Using mt = m^siX^ — X2) and (jA.3p . we can rewrite the first three terms on the right 
hand side of this equation as 

. [2e + m' (A3 - A2)] =~^ = -4^arctanh (^) = -^ In (^—%] • (A.7) 

SA2 SA2 \sX2J 2sA2 \l + il7/ 

Inserting this expression into (|A.6p and making use of (jA.5|) . we obtain as a final result 
(|50p and ([5T|) . where A^ = A2. Real solutions for T exist only when the argument of 
the logarithm is positive, 

which leads to a second inequality to be satisfied by e and m. 

Case (ii) of the above list yields the same result, but with the roles of Ai and A2 
interchanged, i.e. A_l = Ai. Case (iii) can only be satisfied if Ai = A2, but the result 
for T is the same as in (i) or (ii). Case (iv) finally has solutions only for 2e — — m^Aa, 
but the corresponding results for J^ are again the same as in the cases (i) and (ii). 
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